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1 Introduction 



The spin-1/2 antiferromagnetic Heisenberg XXZ chain is one of the most fundamental models 
for one- dimensional quantum magnetism, which is given by the Hamiltonian 

N 

n = J2 (^J^jVi + spl, + ^s^s;^,) , (i.i) 

where 3° = crf/^ with a" being the Pauli matrices acting on the j-th site, A is the anisotropy 
parameter, and N is the number of lattice sites of this system. Here we impose the periodic 
boundary condition Sj+n = Sj. For A > 1, it is called the massive XXZ model where the 
system is gapful. Meanwhile for — 1 < A < 1 case, the system is gapless and called the 
massless XXZ model. Especially we call it XXX model for the isotropic case A = 1. 

The exact eigenvalues and eigenvectors of this model can be obtained by the Bethe Ansatz 
method [Il[2]- We shall give a brief survey of this method below. First we assume the wave 
function |\E') with M down spins in the form 



l<ni<---<nM<N aeSn 



■ M 

i J2 ka(j)nj 



\ni, ■■■,nM), 



A{a) = e{a) JJ (e*('=-(0+'=-0)) + 1 - 2Ae^'=-w) , (1.2) 
Kj 

where a denotes an element of the symmetric group S'„, e(cr) is the sign of the permutation, 
and \ni, um) signifies the state where the spins at the positions rii, ■ ■ ■ ,nM are downward 
and all the other spins directing upward. If the quasi-momenta {ki,k2, - ■ ■ , satisfy the 
Bethe ansatz equations 

^./c,^^TT^_ e +1 2Ae \ (^■ = i,2,...,M), (1.3) 

the wave function (11. 2p becomes the eigenvector of the Hamiltonian (II. ip : 

M 

\ 

4~ ^ 



n\^) = E\m), E = ^ + ^(cosfc^ - A). (1.4) 



In the thermodjTiamic limit N oo, the ground state energy per site cq of the massless 
XXZ model, for example, is computable by analyzing the Bethe ansatz equations 

4 2% smhtcosh(i/t)' ^ ^ 

where we parametrize the anisotropy parameter as A = cos{7tu). Especially in the XXX 
case A = 1, the ground state energy per site (II. 5p can be simplified into 

eo = ^-ln2. (1.6) 
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In the same way other physical quantities in the thermodynamic hmit such as specific heat, 
magnetic susceptibihty, elementary excitations, etc., can be exactly evaluated even at finite 
temperature by the Bethe ansatz method [2]. 

The exact calculation of the correlation functions, however, is still a difficult problem even 
in the simplest case for static correlation functions at zero temperature. The exceptional case 
is A = 0, where the system reduces to a lattice free-fermion model by the Jordan- Wigner 
transformation. In this case, we can calculate arbitrary correlation functions by means of 
Wick's theorem [31 H] . However, there have been rapid developments recently in the exact 
evaluations of correlation functions for A 7^ case also. Below we shall give the historical 
review of them, concentrating mainly on the static correlation functions of the XXX model 
A = 1 at zero temperature for the infinite system ^ 00. Until quite recently only the 
correlation functions within three lattice sites had been known: 

(StS^) = ^ - iCa(l) = -0.1477157268533151 ■ • • , (1.7) 

(•^i^g) = 3^ - ^Ca(l) + Ca(3) = 0.06067976995643530 ■ ■ ■ , (1.8) 

where Ca{s) is the alternating zeta function defined by Ca{s) = V'"-'^! which is 

related to Riemann zeta function ({s) = ^/''^^ Ca{s) = (1 — 2^~'')^(s). Note that 

the alternating zeta function is regular at s = 1 and is given by Ca(l) = In 2. The nearest- 
neighbor correlation function (11.71) is derived directly from the ground state energy (11.61) 
obtained by Hulthen in 1938 |5]. The first nontrivial correlation function (11.81) was derived 
by Takahashi in 1977 via the strong coupling expansion of the ground state energy for the 
half-filled Hubbard model [6j. Note also that another derivation of the second-neighbor 
correlation function (II. 8p was given by Dittrich and Inozemtsev in 1997 [7]. These method, 
however, can not be generalized to calculate further correlation functions on more than four 
lattice sites, unfortunately. 

Meanwhile using the representation theory of the quantum affine algebra Uq{sl2), Ky- 
oto Group (Jimbo, Miki, Miwa, Nakayashiki) derived a multiple integral representation for 
arbitrary correlation functions of the massive XXZ antiferromagnetic chain in 1992 [8, 9j, 
which is before long extended to the XXX case flUi ITT] and the massless XXZ case [T2] . 
More recently the same integral representations were reproduced by Kitanine, Maillet, Ter- 
ras |T3] in the framework of Quantum Inverse Scattering Method. They have also succeeded 
in generalizing the integral representations to the XXZ model with an external magnetic 
field More recently the integral formulas are extended to the XXZ model even at finite 
temperature ! [TU [15]. In this way it has been established now the correlation functions on 
n lattice sites are represented by n-dimensional integrals in general. However, these multiple 
integral representations for correlation functions, though completely exact, have not been 
used widely especially among physicists. It is mainly because we can not evaluate accurate 
numerical values of correlation functions directly from the integral representation. Also it 
had been a puzzle that the exact expression of the correlation function on three lattice sites 
(II. 8p had not been reproduced from the integral representation for a long time. 

The situation changed when Boos and Korepin devised an innovative method to evaluate 
these multiple integrals for XXX chain in 2001 [T6l IT7] . They showed that the integrand 
in the multiple integral formula can be reduced to a canonical form, which allows us to 
implement the integration. This method was at first applied to a special correlation function 
called Emptiness Formation Probability (EFP) [11] which signifies the probability to find a 
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ferromagnetic string of length n: 




(1.9) 



By performing the multiple integrals, the explicit forms of the EFP for up to n = 5 was 
obtained fT6l [IH [18] . 

The Boos-Korepin method was applied to calculate other correlation functions on four 
lattice sites in 2003 [I9]. Especially the third-neighbor correlation function was obtained 
there as 



(^1^4) =^ - 3Ca(l) + yCa(3) - yCa(l)Ca(3) - ^Ca(3)^ - f Ca(5) + yCa(l)Ca(5) 



Other correlation functions on four lattice sites are given in Appendix A. 

Subsequently, the Boos-Korepin method was generalized to XXZ chain both in massless 
and massive regime and all the correlation functions within four lattice sites were obtained 
for general anisotropy [201 [211 [22] . 

In principle, multiple integrals for any correlation functions can be performed by means of 
Boos-Korepin method. However, P(5) is the only correlation function which was calculated 
by this method on five lattice sites, since it is getting critically harder to reduce the integrand 
to canonical form as the integral dimension increases. 

In the course of attempting to obtain further correlation functions, the alternative alge- 
braic method to calculate the EFP was developed by Boos, Korepin and Smirnov in 2003 
[23]. They considered the inhomogeneous XXX model, in which each site carries an inho- 
mogeneous parameter Xj. The homogeneous XXX model corresponds to the case with all 
the inhomogeneous parameters set to be 0. Inhomogeneous correlation functions on n 
lattice sites are considered to be functions of variables Ai, . . . , A^. They derived the func- 
tional relations which the inhomogeneous EFP should satisfy by investigating the quantum 
Knizhnik-Zamolodchikov (qKZ) equations [2H [251 [2S1 [22], the solutions to which are con- 
nected with the inhomogeneous correlation functions [^1 HDl US]- Moreover they suggested 
an ansatz for the form of the inhomogeneous EFP, which consists of only one transcendental 
function u;(A) with rational functions of inhomogeneous parameters Ai, . . . , A^ as coefficients, 
(for the proof of the ansatz and further generalizations, see [281 ISHl [301 [SD [221 [33]). It was 
shown that the functional relations together with the ansatz for the final form completely 
fix the explicit form of the inhomogeneous EFP. In this way the explicit forms of the inho- 
mogeneous EFP for up to = 6 have been obtained, which gives a new result for P(6) in 
the homogeneous limit Aj 0. 

This Boos-Korepin-Smirnov method based on the qKZ equation was generalized to ar- 
bitrary correlation functions in [21J. Actually, in that paper, all the correlation functions 
on five lattice sites are obtained based on the functional relations for the general correlation 
functions (see Appendix A). Especially the fourth- neighbor correlation function is given by 



0.05024862725723524 ■ ■ ■ . 



(1.10) 



h - y^'^^^) + T*^"^^^ " 72Ca(l)Ca(3) - ^Ca(3)^ - 



700 



Ca(5) 




9 



3920 



9 



Ca(l)Ca(7) 



+ 280Ca(3)Ca(7) = 0.03465277698272816 ■ ■ ■ . 



(1.11) 
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The main purpose of this paper is to report further results using the algebraic method in 
[51] . More explicitly we have succeeded in calculating all the correlation functions on six 
lattice sites using this algebraic method. We remark that, if we consider only the two-point 
correlation functions (SiS^), it was already calculated up to the seventh-neighbor correlation 
function (5'fS'|) [35l|36]. The method actually allows us to evaluate some other correlation 
functions such as string correlation functions [37], but not all the correlation functions. 
We also remark that there is a related but different approach for the evaluation of the 
general correlation functions, which is developed recently by Boos, Jimbo, Miwa, Smirnov 
and Takeyama [281 1291 [301 (SB [321 [33] . They have established a compact exponential formula 
for general correlation functions without heavy multiple integrals, which would be also useful 
to evaluate correlation functions analytically. We, however, do not discuss the formula in 
this paper. 

The plan of this paper is as follows. In section 2, the algebraic method to calculate 
general correlation functions is reviewed. In section 3, we present the analytical results for 
some physically interesting correlation functions, such as chiral correlation functions, dimer- 
dimer correlation functions, etc- ■ ■ . In section 4, we show the eigenvalue-distribution of the 
reduced density matrix and calculate the von Neumann entropy. Summary and discussion 
are given in section 5. 



2 Algebraic method to evaluate general correlation func- 
tions based on the qKZ equations 

Below we describe the functional approach to evaluate general correlation functions estab- 
lished in |3l|. Any correlation function can be expressed as a sum of density matrix elements 
Pei,...',lni which are defined by the ground state expectation value of the product of elementary 
matrices: 

pt::i:^{Et''---Ei^'"). (2.1) 

e' e ■ 

where E^^ ^ are 2x2 elementary matrices acting on the j-th site as 

^ hi] ^ hi] 

For example, the density matrix elements on four lattice sites Pl^^lz can be written in terms 
of spin operators as 

P^t^- = (^r + Sl^ St Q - 5|) ^ . (2.2) 

First we rewrite the functional relations for inhomogeneous density matrix elements, which 
follow from the qKZ equations: 
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• Translational invariance 

(Ai + A, . . . , A. + A) = Pt:;::{Xi, ...,Xn), 

• Transposition, Negating and Reverse order relations 

i^:;;;^(Ai,...,A„)=i-^;;;;-(-Ai,...,-A„) 

= p:t:=i{K . . . , A„) = Pt::;:H->^n, ■ ■ ■ , -Ai) 

• Intertwining relation 

y2 -^e'eC! ~ Xj+i)P..,ej+i,e-,... (■ ■ ■ , Aj+i, Xj, . . .) 

=± 



= E 



where R denotes the i?-matrix of the XXX model: 



R{X) 



fRXtW RltiX) RtXiX) RttiX)\ 

RUiX) RXz{X) RUiX) Rtz{X) 

R+XW ^;i(A) Rztix) RzliX) 

\RX+{X) R+Z{X) RZXiX) RZZ{X)J 



(1 








o\ 





A 


1 





A+1 


A+1 





1 


A 





A+1 


A+1 










V 



First recurrent relation 



^::t(A + i,A,A3, 

^'^-•••'^"(A-l,A,A3, 



Second recurrent relation 

■ ,e^{Xi, . . . ,Xj,. 



lim Pe^' 



£1 



■ , A„) 

• , Xn) 



■ ■ , An) 



-Sei-e2^1^2Pt^^'^^'^'X^"^^{X, A3, ... , An) 

-'^4,-4^i4i^:i:er"(A,A3,...,An) 



1 e' e' e' e' 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



• • , Aj_i, Aj+i, . . . , A 



(2.8) 



• Identity relation 

E^'::i(^i'---'^")= E ^:::i(Ai,---,An) 

= -P+,'.'.'.',+ (Al, . . . , An) — P-l'.. '.'- (Al, . . . , An) 

• Reduction relation 

(Al, A2, . . . , An) + PZ±±{X^, A2, . . . , An) = (A2, . . . , An) 



(2.9) 



-1,'+ (''^1' ■ ■ ■ ' An-1, An) + Pell..',e2-i- (Al) ■ ■ ■ ) An-1, An) — Pell-le^-i {Xi, 



, Xn-l) 

(2.10) 
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Additionally it is established that density matrix elements can be written in the form of 



nt:£'(A>,...,A„)=( 



2 

[f] I 

1=1 ct£T„^i j = l 

where T„^/ is a subset of the symmetry group of degree n defined as 

Tn,i = {(re Sn\(y{l) < a(3) < • ■ ■ < a{2l - 1), 

a(2j - 1) < a(2j) for j = 1,2, ■■■,/, 

a{2l + 1) < a(2l + 2) < ■ ■ • < a{n)}. (2.12) 

o"|Ai, . . . , A„) are rational functions of inhomogeneous parameters Ai, . . . , A„ with 
known denominator 

n M nLll-^afe-l - A2fc) n2«+l<fc<i<ri('^fc - Aj) 

JJn^il^Ai, . . . , A„j - ^^j^ — — . [Z.16) 

lll<fc<j<nl^fc - ^j) 

fjAi, . . . , A„) are polynomials of known degree in inhomogeneous parameters Ai, . . . , A„ 
with rational coefficients which have at most the same degree in each variable and also at 
most the same total degree as in the denominator. 

A transcendental function uj{X) in (12.111) is defined by 

^(A) = 2 + = ^ S A''{Ca(2A: - 1) - W^k + 1)}. (2.14) 

fc=l fc=0 

Here note that the identity Ca(— 1) = •~3C(— 1) = 1/4. The following properties of the 
function a;(A) are needed for calculations of the recurrent relations (I2.7l) - (l2.8p : 

o;(zoo) = 0, ^(A±l) = ^p^-^^^^W- (2-15) 

Below we list the explicit forms of all the non-zero elements of inhomogeneous density 
matrices for n = 1, 2: 

p+(AO = Pr(AO = ^, 

(Ai, As) = Prr(Ai, \2) = \ + \uj{\r - A2), 

(Ai, A2) = Pr+(Ai, A2) = i - ^c.(Ai - A2), 

P++ (Ai, A2) = P+1 (Ai, A2) = ^^(Ai - A2). (2.16) 
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It can be easily seen that these satisfy the functional relations (12.3l) - fl2.10l) . By taking 
the homogeneous limit Ai, A2 — > 0, we obtain density matrix elements for the physically 
interesting homogeneous case: 



Pi 



P 



P- 



P++ 



Pi 



I-IUD 
1 

6 



;Ca(l). 



(2.17) 



In the paper [M], it has been argued that the functional relations (I2.3p - fl2.10p together 
with the ansatz for the form of inhomogeneous correlation functions (12. lip completely de- 
termine the unknown rational coefficients in the polynomials Qeii!!!ie^(^, o"|Ai, . . . , A„) and 
therefore the full form of the inhomogeneous correlation functions. As an example, we give 
the explicit form of the density matrix element for n = 3 



g;+l(/ = l,or=(123)|Ai,A2,A3) 
Q;+l(/ = l,a = (231)|Ai,A2,A3) 



g;+l(/ = l,a = (132)|Ai,A2,A3 
2 — A12 



6 " 12^^^' 
1 1 - 

6 + 12^^^' 

11,, 1 , . 

77 + 7:Ai3A23 — t7t(Ai3 + A23J 
00 12 



P 



(Ai, A2, A3) 



12Ai,A 



2 + A23 , 2 + 2A12A32 — (A12 + A32) 

+ —^ ^^23 + H TTTT T t^l3, 



13^23 



I2A21A31 



I2A19A 



(2.18) 



12^32 



where we have used the abbreviations Ujk = uj{XjiS), Xjk = Aj — A^. By taking the homoge- 
neous limit Ai, A2, A3 — i> 0, we obtain the homogeneous density matrix element 



(2.19) 



=^-^Ca(l)+Ca(3). 

In this way all the correlation functions on five lattice sites have been obtained in [M]. In 
this paper we have obtained further results for six lattice sites. There are 24 independent 
correlation functions among them, explicit forms of which are collected in Appendix B. 

Concerning two-point correlation functions (SfS^), on the other hand, an efficient way 
to calculate them has been developed [351 EE] . Namely, without evaluating all the density 
matrix elements, a two-point correlation function can be derived from its generating function 




(2.20) 



through the relation 



1 92 



2dK^ 



JDK 9 JDK, 

n n-1 



JDK. 
-'n-2 



1 

4' 



(2.21) 



Note that the generating function can be considered as a generalization of the Emptiness 
Formation Probability (II. 9p . which is reproduced if we set k = in (I2.20p . From the 
functional relations for general density matrix elements f l2.3p - fl2.10p . it can be shown that 
the generating function P^ satisfies the following closed functional relations: 
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• Translational invariance 

P„^(Ai + A, . . . , A, + A) = P:{Xi, . . . , AO (2.22) 

• Negating relation 

P:{-\u . . . , -A„) = P„«(Ai, (2.23) 

• Symmtry relation 

P„^(Ai, . . . , A„) = P„^(A.(i), . . . , A.(„)), (2.24) 
where a denotes any element of the symmetry group Sn- 

• First recurrent relation 

n-l, An-1 ± 1) — kP^_2{Xi, ■ ■ ■ , Xn-2) 

• Second recurrent relation 

lim P„^(Ai, . . . , A„_i, A„) = i^P^i(Ai, . . . , A„_i) (2.26) 

One can explicitly calculate the generating functions recursively with respect to n from 
these functional relations together with the ansatz for the final form 

I 

, A<^(„)) ]^t^(A^(2j_i) -A^(2j)), 



(2.27) 

where T„^;, Dn,i{Xi, ■ ■ ■ , Xn) and uj{X) are the same as in (12. lip and ^(Ai, . . . , A^) are 
polynomials containing the parameter k. The great advantage of this method is that the 
polynomials Q'^ ^(Ai, . . . , A^) do not depend on the permutation a G Tnj due to the symmetry 

relation (12.241) . unlike the case of general density matrix elements Q'el]Z]l'^{l, a\Xi, . . . , Xn)- 
This fact considerably reduces the amount of bothersome calculation. Actually by using this 
method, up to the seventh- neighbor correlation function {SfS^) have been calculated [36] . 
Calculating all the density matrix elements is much harder. We, therefore, have succeeded 
only up to six sites yet. 

3 Applications to physically interesting correlation func- 
tions 

In this section, we discuss several physically interesting correlation functions, which can be 
evaluated exactly from our results of all the density matrix elements for n = 6. 
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2.25 



P^(Ai, . . . , A„) 



K 



1=1 cr£Tn,i 



A'^l{Xi, . . . ,Xn) 



(1)' 

■ ■ ■ y^n) 

Dn,l{Xl, . . . , Xn) ' 



3.1 Chiral correlation function 

First let us consider the vector chiral correlation function 

((^i X S2) ■ [S, X = 6 {{S^.S^S^S]^,) - , (3.1) 

which measures a chirality of the spin alignment. It has been observed that they have simple 
expressions [311 EH] compared with the other correlation functions on the same lattice sites 

(Si X ^2) ■ (5^3 X 54)) = i (Ca(l) - Ca(3)) = -0.1041977484048752- ■ ■ , (3.2) 

(5i X S2) ■ (^4 X 55) ) = i (Ca(l) - Ca(3)) " ^ (Ca(3) - Ca(5)) = "0.01597638205835821 ■ ■ 

(3.3) 

Here we have newly obtained the vector chiral correlation function for six lattice sites, which 
has also a simple form though it contains some quadratic terms 

((^1 X S2) ■ xSe)) = ^Ca(l) - yCa(3) + 9Ca(5) - yCa(7) + ^Ca(l)Ca(3) 

4 39 4 98 

+ 3Ca(3)^ - yCa(l)Ca(5) - -C(3)Ca(5) + yCa(l)Ca(7) 

= -0.01774606473688137 ■ ■ • . (3.4) 

It is interesting to note that the numerical values of these are all negative and their absolute 
values are oscillating. 

Furthermore we obtain the scalar chiral correlation function defined as below, which 
reveals an intriguing factorization: 

([(51X52) -53] [(54X5^5)-56]) 

at I C'J: QV QZ qx qV qz\ , / QX qV qz ny qz qx\ , / nx ny nz qz qx qS/\ 

— Ul^\0]^ 0203040505/ -t- \0]^ 0203040505 / -t- \0;^ 0203040505/ 

/ QX oy QZ QZ Qy QX\ r\ I QX QV QZ QX QZ QVX \ 

\0;^ 0503040505/ ^^0^^0203040505/; 

= liUl) - l/4)(Ca(7) - Ca(5)) - ^(Ca(3) - Ca(l))(Ca(5) - Ca(3)) 

= 0.008133862120680087 ■ ■ • . (3.5) 

Note that the numerical value in this case is positive. These numerical values may indicate 
a classical picture of the spin alignment for the antiferromagnetic ground state of quantum 
Heisenberg chain. 



3.2 One-particle Green function 

By the Jordan- Wigner transformation 

i-i j-i 
5+ = Ha - 24c,)c„ Sr = l[{l-2clc,)cl, (3.6) 

k=l k=l 

the XXX model is transformed into the isotropic spinless fermion model. Let us consider 
the one-particle Green function {c\cn)f for this model. Here the bracket (■■■)/ means the 
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expectation value in the half-filled state of the spinless fermion model. The first-neighbor 
one-particle Green function (c{c2)/ is directly obtained from the Hulthen's result (11.71) as 

(c|c2)/ = {S^S^) = 2{S^S^) = ^ - ^Ca(l) = -0.2954314537066302 ■ • ■ . (3.7) 
The first non-trivial result for this one-particle Green function has been obtained in [19] 

{clc,)f = MStS^S^S^) = 8{stsis',s:) 

= ^ - 2Ca(l) + ^Ca(3) - yCa(l)Ca(3) - |Ca(3)^ - ^Ca(5) + yCa(l)Ca(5) 

= 0.08228771668643698 • ■ • . (3.8) 

because (c|c„)/ = if n is odd. Here we have newly obtained the fifth-neighbor one-particle 
Green function 

{clc,)j = l6{S^S^,SIStSIS^) = 32{S'tS^,SISlSISI) 

1 10, , , „ . . X 3608, . X. . X 90832, , ,0 100288, , ,0 

= yCa(l) + 74Ca(3) - ^Ca(l)Ca(3) - ^^Ca(3)' - ^^Ca(3)^ 

- ^a5) + l^ai)a5) + ^a3)a5) + ^ai)a3)a5) 

- '-^cM^) - ^as)^ + ^c(i)a5)^ + ^a3)a5)^ 

490880,,^,, 89918,,,, 308392,,,,,,,, 623128, ,„ 
+ + ^g-^"*^) —CaWUT) + — ^Ca(3)C(7) 

_ lH|H5c.(i)c(3)c(7) - '-^um^m + ^c.(5)c.(T) 
_ 2Hi5a3)U5)c.(7) - ?^u7r . '-^um^f - ^cp) 

128928, , , , 196304, , , , 645568, ,,,,,,,, 

+ ^^Ca I Ca 9 3 Ca 9 + _— 1 (a 3 (a 9 

5 3 15 

98176, , ,o, , , 135104, , , , 196352, ,,,,,,,, 

+ Ca(3)'Ca(9) + Ca(5)Ca(9) ^Ca(l)Ca(5)Ca(9) 

= -0.04497471675792834 ■ ■ ■ . (3.9) 

It will be an interesting problem to study the asymptotic behavior of one-particle Green 
function in detail. Conformal field theory predict that 

{c\c^)f ~ cos {kFx) a;"^/^ kp = n/2, (3.10) 



omitting the logarithmic correction [52]. Unfortunately our new exact results are not suffi- 
cient to confirm the asymptotic formula. 

3.3 Dimer-dimer correlation function 

Next let us consider the dimer-dimer correlation function defined as 

Dn = (^("S^i ■ ^'ij (^^n ■ Sn+i^^ = 3{S^S2S^S^^-^} + 6{S^ S2 S^S^_^_^) , (3.11) 
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which should asymptotically coincide with the square of the ground state energy 



n— »oo 



1 



lim D„ = = ( - - Ca(l) ) = 0.1963794236382287 ■■■ . (3.12) 



From our results of all the density matrix elements for tt, = 6, we can exactly evaluate the 
dimer-dimer correlation functions Dn for up to n = 5: 

e; 4 ^10 

D,-el= -C(3) - -Ca(l)Ca(3) - Caiif - -^(5) + yCa(l)Ca(5) - U^f 

= 0.06082478294036410 ■ ■ ■ , (3.13) 

D,-el= ^Ca(3) - yCa(l)Ca(3) - yCa(3)^ - yCa(5) + ^Ca(l)Ca(5) - 5^(3)^(5) 

- f Ca(5)^ + ^Ca(7) - ^Ca(l)Ca(7) + 35^(3)^(7) - C(l)' 

= -0.02773785800119889 ■ ■ • , (3.14) 

D,-el = ^C(3) - 10C(l)Ca(3) - ^-fCaiSf - ^-fCaiSf - |^C(5) + ^C(l)Ca(5) 

196 640 904 

- — C.(3)C(5) + — C(1)C(3)C,(5) - ^Ca(3)'C.(5) - 1381C.(5)^ 

+ ^ai)c.(5)= + '-fusKinr + ^a5)' + 'fur) - '^cmcp) 

3 3 9 6 3 

+ 3038C(3)C(7) - ^Ca(l)Ca(3)Ca(7) - l^Ca(3)^Ca(7) + ^Ca(5)Ca(7) 

- ^Ca(3)Ca(5)Ca(7) - ^Ca(7)^ + ^Ca(l)Ca(7)^ - ^Ca(9) + 798C.(l)Ca(9) 

1 4224 4368 

- ^Ca(3)Ca(9) + 1904Ca(l)Ca(3)Ca(9) + -^Cai3)%i9) + 2156Ca(5)Ca(9) 

5 5 

- 2912Ca(l)Ca(5)Ca(9) - Ca(l)' 

= 0.01892813120084483 ■ ■ ■ , (3.15) 
The difference D„ — Cq also should decay algebraicly. 



3.4 AP{n): a probability to find an antiferromagnetic string of 
length n 

Emptiness Formation Probability P{n) (11.91) represents the probability to find a ferromag- 
netic string of length n. Similarly we may consider a probability to find an antiferromagnetic 
string of length n, which we denote by AP{n) 

AP{n) = p_^r^z: + Pr+r+;:. (3.16) 

From our results we can exactly evaluate the AP[n) for up to n = 6 

AP{2) = P+- + Pzt = \ + ^Ca(l) = 0.7954314537066302 ■ • • , (3.17) 
^P(3) = P+I+ + Pzt- = ]:- ^Ca(l) + Ca(3) = 0.6061112236630655 ■ ■ ■ , (3.18) 
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AP{A) = p+r+r + Pitit 

= ^ - ^Ca(3) + ^Ca(l)Ca(3) + ^Ca(3)^ + ^Ca(5) - ^Ca(l)Ca(5) 

= 0.4938083479102196 ■ ■ ■ , (3.19) 

AP{5) = p+r+r+ + Pitt 

- f Ca(3)Ca(5) - ^Ca(5)^ + ^Ca(7) - ^Ca(l)Ca(7) + ^Ca(3)Ca(7) 

= 0.3943245947356898 ■ ■ ■ , (3.20) 

AP{6) = Pttt + Pttt 

= 4 - iCa(l) + |Ca(3) + f Ca(l)Ca(3) + ^Ca(3)^ 

+ ^^3)3 + ^a5) - ^ai)a5) - f a3)a5) 

- ^C.(l)C.(3)C.(5) + ^C.(3)^C«(5) + ^a(5)^ - ^C.(1)C.(5)^ 

- ^a3)a5)^ - i^a5)3 - ig^ar) + ^uDur) 

9576 26152 ^ 1868 

^Ca(3)Ca(7) + ^-Ca(l)Ca(3)Ca(7) + ^^Ca(3) Ca(7) ^Ca(5)Ca(7) 

+ ^a3)U5)a7) + i^ar)^ - '-^cMU7y- . f us) 

- ^fiui)c.(9) + ^a3)U9) - ^ai)a3)U9) 

19868, , , X 27028, . x> . x 39736, , , , , 

- ^^Ca(3)^Ca(9) - ^^Ca(5)C(9) + ^^Ca(l)Ca(5)Ca(9) 

= 0.3239037769698205 • • • . (3.21) 



It is known that the EFP P{n) shows a Gaussian decay jlQl HI] 

P{n) ~ An-^C""'. (3.22) 

On the other hand, the analytical asymptotic behavior of AP{n) has not been stud- 
ied. However, Y. Nishiyama noticed that it decays roug hly AP{n) ^ [AP{2)f from his 
numerical data before [12] . We compare the formula with our exact data in Table [H 



4 Reduced density matrix and entanglement entropy 

At zero temperature, the density matrix for the infinite system can be written as 

PT = |GS)(GS|, (4.1) 

where |GS) denotes the antiferromagnetic ground state. Let us consider a finite sub-chain of 
length n in the infinite chain, the rest of which can be regarded as the environment (Figure 
[1]). The reduced density matrix for this sub-chain is obtained by tracing out the environment 
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Table 1: Asymptotic behavior of AP{n) 



n 


AP{n) Exact 


r * -1—1 / ^ \ tTJ — - 1 

[AP{2)] 


2 


0.795431 


0.795431 


3 


0.606111 


0.632711 


4 


0.493808 


0.503278 


5 


0.394325 


0.400323 


6 


0.323904 


0.318430 



Environment 



System 



OHDH>CHDH> - - - <h:><><><> - - 



12 3 



n-1 n 



Environment 



Figure 1: Finite sub-chain of length n in the infinite chain 



from the infinite chain: 

pn = tr£;PT 



(4.2) 



From our results, we have computed all the eigenvalues cu^ (a = 1, 2, ■ ■ ■ , 2") of the reduced 
density matrix p„ for up to n = 6, which are shown in Figure [2l 

We have found that the smallest eigenvalue a;2" is, for any n, the emptiness formation 
probability Pin). This is in + l)-fold degenerate since the reduced density matrix has the 
block-diagonalized form by the magnetization of the sub-chain = 'YTj=i^'j ~ Ylj=i^j^ 
each of which has the non-degenerate smallest eigenvalue P{n). Or in other words, it is a 
consequence of SU{2) invariance of the Hamiltonian. 

From these results we can calculate the von Neumann entropy (Entanglement entropy). 

2" 

S{n) = -trp„ log2 Pn = - ^ UJa log2 (^a, (4.3) 

0=1 

which is considered to measure how the ground state is entangled [13] . The exact numerical 
values of S{n) up to n = 6 are shown in Table [21 It can be seen that as n grows the von 
Neumann entropy is increasing in small steps. The asymptotic behavior of this entropy is 
discussed in |13]. In the massive region A > 1, the von Neumann entropy will be saturated 
as n grows, which means the ground state is well approximated by a subsystem of a finite 
length corresponding to the large eigenvalues of reduced density matrix. On the other hand, 
in the massless case — 1 < A < 1, the conformal field theory predict that the von Neumann 
entropy shows a logarithmic divergence [l4j, of which the explicit form for the XXX case 
A = 1 reads 

S{n) ~ ^ log2n + k. (4.4) 
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10-^ 



10"^ 



10 ^ 



10-^ 



10 ^ 



10"^ 



SB 


1 1 1 






n=l 


■ 


▲ A 






• • 


n=3 


o 




▲ ▲ ▲ ▲ 






oooo n=5 


A 




▲ ▲ 






▲ ▲ ▲ ▲ 






A A 






A A A A 






A A 






A A A A 






A A A A 

1 1 1 


A A 



10 20 30 

a 

Figure 2: Eigenvalue-distribution of density matrices 
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Table 2: von Neumann entropy S{n) of a finite sub-cliain of length n 



S{1) 


S{2) 


^(3) 


^(4) 


1 


1.37585732628 


87466 1.5824933209573855 


1.7247050949099274 



S{5) S{6) 
1.833704916848315 1.922358833819333 



We estimate the numerical value of the constant term k = 1.06209 by the extrapolation 
S{n) — I log2 n = Co + Ci/n + 02/^^. Our exact results for up to n = 6 agree quite well with 
the asymptotic formula as shown in Figure [31 

5 Summary and discussion 

We have succeeded in obtaining the analytical forms of all the density matrix elements on 
six lattice sites for Heisenberg chain using the algebraic method based on the qKZ equations 
developed in [3l]. We have obtained several physically interesting correlation functions, 
such as chiral correlation functions, dimer-dimer correlation functions, etc- ■ ■ . Subsequently 
we have calculated all the eigenvalues of the reduced density matrix p„ for up to n = 6. 
We observe that the smallest eigenvalue is the emptiness formation probability P{n) and 
is (n + l)-fold degenerate. Of course, it may be more desirable if we could describe the 
eigenvector corresponding to the largest eigenvalue. From these results we have computed 
the von Neumann entropy, which shows a good agreement with the asymptotic formula 
derived via the conformal field theory. 
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Appendix A Known results for up to n = 5 

In this appendix we present all the independent correlation functions for up to n = 5 except 
for the two-point correlation functions (SfS^) given in Introduction. First we give the 
correlation functions on four lattice sites [19] : 

1 1 qi 2 1 7 2 

{S^^S^.S^S!) =— + -Ca(l) - Y^Ca(3) + -Ca(l)Ca(3) + -Ca(3)^ + -Ca(5) - -Ca(l)Ca(5) 

=0.02750969925180030 ■ ■ ■ , (A.5) 
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n 



Figure 3: von Neumann entropy S{n) of a finite sub-chain of length n 



{SIS^SISD ~ - lCa(l) + ^Ca(3) - ^Ca(l)Ca(3) - ^Ca(3)^ - ^Ca(5) + ^C(l)Ca(5) 



- 0.007080326815007911 



(A.6) 

{SIS^SISD ~ - iCa(l) + ^Ca(3) - ^Ca(l)Ca(3) - ^C(3)^ - ^Ca(5) + ^Ca(l)Ca(5) 

=0.01028596458580462 ■ ■ • , 

/ QZ QZ QZ QZ\ / QX QX QZ QZ\ i / QX QZ QX QZ\ i / QX QZ QZ QX\ 

\>-'l>J2'-'3'-'4/ — \'-'l'-'2'-'3'-'4/ "T X*-"! '-'2'-'3'-'4/ "r Wl '-'2 '-'S '-'4 / • 



(A.7) 
(A.8) 



Note that on the antiferromagnetic ground state without magnetic field, the correlation 
functions with an odd number of spin operators vanish. Considering the isotropy of the 
Hamiltonian (11.11) . we can see that the independent correlation functions on four lattice sites 
are restricted to the above ones and the two-point correlation function (SfS^). 

Correspondingly the independent correlation functions on five lattice sites are [31] 



1 l./x 517^,, 25^,,^,, 203^,,. 305,,, 
172 4 40 35 49 

- — Ca(l)Ca(5) + -Ca(3)C(5) + yCa(5)^ - ^Ca(7) + yCa(l)Ca(7) 

- — Ca(3)Ca(7) = -0.009892435084700294 • • • , 

3 



(A.9) 
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{S^S^2SISI) ~ lCa(l) + ^Ca(3) - 6C(l)Ca(3) - ^Ca(3)^ - yCa(5) 

+ ^Ca(l)Ca(5) - gCa(3)Ca(5) " ^Ca(5)^ + ^Ca{7) - 42C„(l)Ca(7) 

77 

+ ^Ca(3)Ca(7) = 0.002788973399542967 • • • , (A.IO) 

(^r^l^a^^S^) =^ - ^Ca(l) + ^Ca(3) - f Ca(l)Ca(3) - ^Ca(3)^ - ^Ca(5) 

+ ^Ca(l)Ca(5) - gCa(3)Ca(5) - ^Cai^ f + f Ca(7) - 42C(l)Ca(7) 

+ ^Ca(3)Co(7) = -0.005056660894819286 • • • , (A.ll) 

/ nz QZ QZ\ / QX QX QZ QZ\ _, / QX QZ QX QZ\ _, / QX QZ QZ QX\ A 1 

\*-'l^2^3^5/ ~Wl^2^3^5/ I Wl^2^3^5/ ' Wl^2^3^5/' 



11 41 Q 70 221 

=— + -C„(l) - — Ca(3) + -Ca(l)Ca(3) + — Ca(3)^ + ^Ca(5) 

- ^Ca(l)Ca(5) + yCa(3)Ca(5) + ^Ca(5)^ - ^Ca(7) + -^Ca(l)Ca(7) 

91 

- ^Ca(3)Ca(7) - 0.01857662093837095 • • • , (A.13) 

{SIS^SIS^ ~ - ^Ca(l) + ^Ca(3) - ^Ca(l)Ca(3) - ^Ca(3)^ " ^Ca(5) 

+ ^Ca(l)Ca(5) - y Ca(3)Ca(5) " ^Ca(5)^ + ^Ca(7) - ^Ca(l)Ca(7) 

154 

+ ^Ca(3)Ca(7) = -0.001089368972914637 ■ ■ ■ , (A.14) 
1 1 4Q7 119 1 1 09 77 

{SISISISD - -C(l) + ^Ca(3) - — Ca(l)Ca(3) - ^Ca(3)^ - yCa(5) 

+ ^Ca(l)Ca(5) - yCa(3)Ca(5) - ^Ca(5)^ + |^Ca(7) - ^Ca(l)Ca(7) 

154 

+ Ca(3)Ca(7) = 0.001573361370145065 • • • , (A.15) 

3 

{SIS'^SISD ={SISISISI) + {SIS^SISD + (S^^l^l^g^). (A.16) 



Appendix B Explicit form for all the correlation func- 
tions on 6 lattice sites 

In this appendix the exphcit results for all the correlation functions on 6 lattice sites are 
shown. There are 24 independent correlation functions among them. Any other correlation 
functions can be computed by the following relations: 



{sisisisD = {s^.sisisD + {s^.sisisD + {s^.sisisi), (b.i) 

/ QZ QZ QZ QZ\ / QX QX QZ QZ\ , / QX QZ QX QZ\ , / QX QZ QZ QX\ f-Ti r)\ 

Wl^2^5^6/ — Wl*-'2*-'5*^6/ ^2^5*^6/ ^ \^1^2^5^6/' {^•^) 

I QZ QZ QZ QZ\ / QX QX QZ QZ\ , / QX QZ QX QZ\ , / QX QZ QZ QX\ /"p Q\ 

{StSISIS^,) = ^1^6 ) + {S^SISIS^,) + {S^.SISIS^,) , (B.4) 
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/ QX qV qV qz qz\ 


/ QX QX QX QX QV Qy\ 

— \'-'l'-'2'-'3'-'4'-'5'-'6/ 


/ QX QX qV qz qV qz\ 
~ X*-"! '-'2'-'3'-'4'-'5'-'6/ 


/ QX QX qV qz qz Qy\ 
~ \'-'l'-'2'-'3'-'4'-'5'-'6/' 


(B.5) 


1 QX CU QX QZ QV QZ\ 


/ QX QX QX Qy QX Qy\ 

— \0]^ O2 O3 O4O5 Oq/ 


/ QX- QX Qy QZ Qy QZ\ 

\Oi 0203040505/ 


/ QX Qy QX QZ QZ Qy\ 

\0^ O2O3 040505/, 


(B.6) 


/ QX Qy QP qz qz QX\ 
D2 Og D4 Dg / 


/ QX qx QX qy Oj/ qX\ 

\0-|^ O2 O3 O4 Oc^Oq / 


/ OX- qx qy qz qZ qy\ 

^2 3 4 5 6 / 


/ qx qy qX qZ qz qy\ 
\0-|^ O2 O3 O4O5 O5 / , 




/ QX QV QZ QX QV QZ\ 

'-'2*-'3*-'4*-'5*-'6/ 


/ QX QX Qy QX QX Qy\ 

— \'-'l'-'2*-'3*-'4*-'5*-'6/ 


/ QX QX QV QZ QZ Qy\ 

'-'2*-'3*-'4*-'5*-'6/ 


/ QX Qy QZ QX QZ Qy\ 
\OlO2O3O4O5Og;, 


(B.8) 


/ C^; qV QZ qD qz qx\ 

0203040500/ 


/ QX QX QV QX QV QX\ 


/ QX QX qV qz qV qz\ 
— \Di 0203040505/ 


/ qv qz QX QZ Qy\ 
— \0]^ 0203040505/, 


(B.9) 


/ QX QtJ QZ QZ QX Qy\ 


1 QX Qy QX QX QX Qy\ 

— \*-'l*-'2^3^4'-'5*-'6/ 


/ QX Qy QX QZ QZ Qy\ 
\Ol O2O3 O4O5OQ/ 


/ QX Qy QZ QX QZ Qy\ 
— (^0^^0203040505/, 


(B.IO) 



/ QX QX QV QV QX QX\ / QX Qy QZ QZ Qy QX\ , / QX QX QX QX QV Qy\ 1 / QX Qj/ QX QX QX Qj/\ 

^0^^0203040505/ — (^0^^0203040505/ \0;^02 03 04 05 05/ -|- \0;^ O2 O3 O4 O5 O5 / 

/ QX QX Qy QZ Qy QZ\ / QX QX Qj/ QZ QZ Qy\ I QX Qy QX QZ QZ Qy\ I QX Qy QZ QX QZ Q'i/\ ('R 1 1 X 

\Ol O2O3O4O5O5/ ^0^^0203040505/ — ^0^^0203040505/ — \0;^020304 0505/, 

/ QX Qy QX QX Qy QX\ / QX Qy QZ QZ Qy QX\ , / QX QX QX Qj/ QX Q|/\ 1 / QX QX Qj/ QX QX q3/\ 

^0^^0203040505/ — ^0^^0203040505/ i- \0^020304 0505/ -|- \0;^ O2 O3 O4 O5 O5 / 

/ QX QX Qj/ QZ Qy QZ\ / QX QX Q?/ QZ QZ Q'f/\ / QX Qj/ QX QZ QZ Q'f/\ / QX Q'f/ QZ QX QZ QS/\ /-p 1 0\ 

yO]^ O2 O3 O4 O5 O5 / y O-j^ O2 O3 O4 O 5 O5 / ^0-1^0203040505/ yO-|^ O2 O3 O4 O5 O5 / , ySD.XZiJ 

I QX Qy Qy Qy Qy QX\ / qx Qy qz qz Qy QX\ , / qx QX QX Qy Qy QX\ I / QX QX Qy QX Qy QX\ 

^0^^0203040505/ — \0;^ O2O3O4O5O5/ -|- \0]^ O2 O3 O4O5O5 / -|~ ^0^^0203040505/ 

/ QX QX QV QZ Qy QZ\ / QX QX QV Q Z Q Z Qy \ / QX Qy QX Q Z Q Z Qy \ / QX Qy Q Z QX Q Z Qy \ fT3TQ\ 

— \O1O2O3O4O5OQ; — \O1O2O3O4O5Og; — \O1O2O3O4O5Og; — \O1O2O3O4O5Og;, {D.IO) 



I QX QX QX QX QX QX\ / QX QV QZ QZ QV QX\ _i_ / QX QX QX QX Qy C2/\ _|_ / QX QX QX QV QX Of \ 

\O1O2O3O4O5O5; — \0]^ 0303040505/ -+- \0;^ 0303040505/ -+- ^0^0203040505/ 

I / QX QX QX Qy Qy QX\ 1 / QX QX Qy QX QX Q2/\ _|_ / QX QX QV QX Qy QX\ 1 / QX Qy QX QX QX QV 

"T \Oi0203040505/ -|- \0;^ O2 O3 O4 O5 O5 / -|- ^0^^0203040505/ -|- \0;^ O2 O3 O4 O5 O5 

/ QX QX Qy QZ Qy QZ\ / QX QX QV QZ QZ Q2/\ / QX Q!/ QX QZ QZ Q2/\ / QX Q2/ QZ QX QZ Of \ 

— \Oi O2O3O4O5O5; — \Oi 0203040505/ — \O1O2O3O4O5Og; — \0]^ 0303040505/ 

Below we shall give the explicit results for the independent ones. 



;B.14) 



1 25 800,,^, 1192^,^,^,^, 15368^,^,2 



-Ca(3) - ^Ca(l)Ca(3) - 
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-Ca(3)^ 



^a5) + ^ai)c(5) - f c(3)a5) 

+ 3648C„(l)Ca(3)Ca(5) - ^Ca(3)'Ca(5) - ^Ca(5)' 



608Ca(3)= 



^5!|^C.(1)C(5)^ + ^C.(3)C.(5)= 



3 . 3 

+ 56952C„(3)C„(7) - ^^Ca(l)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 



35392 

3 

66640 



Ca(3)\„(7) + 7840C„(5)C„(7) - 8960C„(3)C„(5)C„(7) 
C„(7)2 + 31360Ca(l)Ca(7)' - 686C„(9) 



+ 18368C„(l)Ca(9) - 53312Ca(3)Ca(9) + 35392Ca(l)Ca(3)Ca(9) 

+ 16128Ca(3)\a(9) + 38080Ca(5)Ca(9) - 53760Ca(l)Ca(5)Ca(9) 

= -0.03089036664760932 • • • 



(B.15) 
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1 1 2^Q dfi"^ IQ^l 

{SfSJSJSD = — + -C.(l) - -^U3) + — C(1)C.(3) + -j^CM' 

2404, 1348,,., 8918,,,,,,,, 3127,, „,,,., 

+ + - ^a(l)C.(5) - — C.(3)C.(5) 

- 25|5ai)U5)^ - '-fuw - ^c(5)» - f c(T) 
+ 5^a(i)c.(T) - ?f c.(3)c.(7) + 5^a(i)c.(3)aT) 

+ ifia3)^av) - ^c.(5)a7) + ^a3)U5)c.(7) 

- ^aDcw + H2!lla(3)a(9) - ^ai)c.(3)a(9) 

12656, ,^,2. 16576,,^,,,^, 25312,,,,,,^,,,^, 

- ^^Ca(3)\a(9) - ^^Ca(5)Ca(9) + ^^Ca(l)Ca(5)Ca(9) 

= 0.003681507672875026 • • • (B.16) 



- - ^^"^'^ + ^Ca(l)Ca(5) + |^C(3)Ca(5) 

12512 6292 85687 , 

+ ^^Ca(l)Ca(3)Ca(5) - ^^Ca(3) Ca(5) ^^«(^) 

25168, ,,,, ,^,2 2864, ,^,2 4352,, 1771,, 

+ -^Ca{l)U5Y + -g-Ca(3)Ca(5)' + -^Cai^f + ^Ca(7) 

89929,,^,^,^, 191513,,^,,,^, 44044,,^,^,^,^,^, 

- ^^Ca(l)Ca(7) + ^^Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

40096, ,^,2. 26866, 30464 , , , 

- ^^Ca(3)'Ca(7) + ^^Ca(5)Ca(7) - ^^Ca(3)Ca(5)Ca(7) 

_Z^c„(rr^i^ai)a7)-^c^^^ 

+ ^c(i)a9) - ^c(3)a9) + ^ai)a3)a9) 

30464, , , , 41944, , , , 60928, . ^> / ^. . ^ 

+ ^^Ca(3)\„(9) + ^^Ca(5)Ca(9) - ^^Ca(l)Ca(5)Ca(9) 

= -0.001116347734065082 • • • (B.17) 
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2452,, 10201,,., 9524,,,,,,,, 37,,„,,,., 

- ^C«(3) - + ^C.(1)C(5) + ^C.(3)C(5) 

+ i^a(i)c„(3)c„(5) - li5c.(3)^c.(5) - '-^u^f 

- '-^uw) + ^a3)c.(7) - ^ai)a3)a(7) 



^^^w^N2 42728 _,,__2 3507 

Ca{7Y + — — Ca(l)Ca(7)' - ^Ca(9) 



5 15 50 

+ ^C(l)Ca(9) - i^Ca(3)Ca(9) + ^G(l)Ca(3)Ca(9) 

36624^,^,2,.,^, 16968,,,,,^, 24416 

+ ^^Ca(3)\a(9) + ^— C„(5)Ca(9) - ^— C„(l)C„(5)C„(9) 

= 0.003069653070471227 • • • (B.18) 



(^r^i^i^i) = 4 + ^c(i) - ^c(3) + ^c(i)Ca(3) + ^c(3)^ 

- ^c(i)a3)c(5) + ^a3)^a5) + ^a5)^ 

215264 2728 -^^^^^^2 4144 3 61327 

271264 _ 61474 376712 _ 

+ -^Ca{l)Ca{7) - Ca(3)Ca(7) + ^^Ca(l)Ca(3)Ca(7) 

38192^ ,^,2^ 79436, 29008 , , , 

+ ^^Ca(3)'Ca(7) - ^^Ca(5)Ca(7) + ^^Ca(3)Ca(5)Ca(7) 

23373, ,„,2 101528, , ,„,2 11977, , , 

- ^C(1)C(9) + ^C(3)Ca(9) - ^C(1)C(3)C(9) 
87024, , ,0, , , 40068, , , , 58016, ,,,,,,,, 

- ^^Ca(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)Ca(9) 

= 0.02384723373803033 • • • (B.19) 



21 



27224, 38389, 35666, ,,,, 18902, 

41491K,^__ 284228, 580888,,^,,,^,,,^, 

- ^^Ca(l)Ca(7) + ^^Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

176624 ^ , 122164 ^ 134176 ^ 

4^^'^^^) + 45 Ca(5)Ca(7) ^^Ca(3)Ca(5)Ca(7) 

324968, _ 2 469616, _ 2 18123,,^, 

- -^Cai7r + ^^Ca(l)Ca(7)^ - ^Ca(9) 

34748, 1343944, 176624, 

+ Ca(l)Ca(9) ^Ca(3)Ca(9) + ^^Ca(l)Ca(3)Ca(9) 

134176, ,^,2. 185696,,^,,,^, 268352,,,,,,^,,,^, 

+ ^^Ca(3)\„(9) + —^Ca{5)Ca{9) - ^^Ca(l)Ca(5)Ca(9) 

= -0.001365096861940014 • • • (B.20) 



{S^SISISD = ^ - ^Ca(l) + ^Ca(3) - ^Ca(l)Ca(3) " ^C(3)^ 

- ^^«^^)^ - ^ ^^(1)^(5) + ^^(3)^(5) 

+ ^ai)a3)a5) - ^a3)^c(5) - ^(.(5^ 

331936, ,,,, ,^,2 12616, ,^,2 19168,,^,, 30961 , 

+ ^^Ca(l)Ca(5)' + -g-Ca(3)Ca(5)' + -^Cai^f + -^Ca{7) 

46109,,^,,,^, 284228,,^,,,^, 580888,,^,,,^,,,^, 

- Ca(l)Ca(7) + ^^Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

176624, ,,2, 122164,,,, _ 134176 , ,,,, _ 

- ^^Ca(3)'Ca(7) + ^^Ca(5)Ca(7) - ^^Ca(3)Ca(5)Ca(7) 

-^a7)^ + ^c(i)a7)-i^c(9) 

+ ^c(i)a9) - l^c(3)a9) + ^c(i)a3)a9) 

134176, , ,o, , , 185696, , , , 268352, ,,,,,,,, 

+ ^^Ca(3)\a(9) + ^^Ca(5)Ca(9) - ^^Ca(l)Ca(5)Ca(9) 

= 0.001592580594206881 • • • (B.21) 
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+ -J^U^) + - -^Ca(l)Ca(5) - ^Ca(3)Ca(5) 

5408^ 8272^,^,2^,, 39289^,^,2 

- ^Ca(l)Ca(3)Ca(5) + ^Ca(3)'Ca(5) + ^g-CalS)' 

- ^c(i)c(5)^ - if^c(3)a5)^ - l^a5)3 - i|^a7) 
+ ^ai)a7) - ^a3)a7) + ^^1)^3)^7) 

+ ^C.(3)\.(7) - ^Ca(5)C.(7) + ^Ca(3)Ca(5)C„(7) 
33124^,^,2 46648 _^__2 1729 

- 714C(1)C(9) + ^^Ca(3)Ca(9) - ^Ca(l)Ca(3)Ca(9) 

13328, ,^,2,^ /^x 18928,, 26656 

- ^^Ca(3)\a(9) - ^^Ca(5)Ca(9) + ^^Ca(l)Ca(5)Ca(9) 

= 0.009188091173609528 • • • (B.22) 



1 11 254 144 6517 

(SfS|SfS|) = — - -C.(l) + — C<.(3) - — C<.(1)C.(3) - -J^CM' 

_ I^a^? - ^a5) + 5|!c.(i)c.(5) . f c.(3)c.(5) 

91712, ,,,, ,^,2 3484, ,^,2 5296, ,^,0 9793, 
106547,,^,,,^, 47236,,^,,,^, 160496,,,,,,^,^,^, 

- ^^Ca(l)Ca(7) + -^Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

- ^C«(3)^C«(7) + ^C«(5)C«(7) - ^C«(3)C.(5)Ca(7) 
91826, ,„,2 129752, , ,„,2 1848, , , 

+ 1778Ca(l)Ca(9) - ^^C(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 

37072 52472 74144 

+ ^Ca(3)\„(9) + ^Ca(5)Ca(9) " ^Ca(l)Ca(5)Ca(9) 

= -0.003158274321296133 • • • (B.23) 
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(5f5;s|5S = ^ - |Ca(l) + ^CO) - ^C.(l)Ca(3) - ^U3f 

- ^-^Uif - ^C,.(5) + ^C(l)Ca(5) - J^C.(3)C.(5) 

+ ^C„(l)C.(3)a5) - If C„(3)^a5) - ^C.(5)^ 

+ ^U1)C(5)^ + 5|lc.(3)a5)^ + 5|5c„,5)= ^ ?^U7) 
108892^ 4755U 160496^ ,^ _ _ 

- ^^Ca(l)Ca(7) + -g-Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

48776 , 3682 37072 

4^^«(2) + ^-Ca(5)Ca(7) ^Ca(3)Ca(5)Ca(7) 

91826^,^,2 129752 _,^_2 1848 

+ 1778Ca(l)Ca(9) - ^^Ca(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 

75 15 
37072^,^,2,,^, 52472,, 74144,,,,,,,,,^, 

+ ^^Ca(3)\a(9) + ^^Ca(5)Ca(9) - ^^Ca(l)Ca(5)Ca(9) 

= 0.005680615367538651 • • • (B.24) 



{S^.SISISD = ^ + ^Ca(l) - ^Ca(3) + ^Ca(l)Ca(3) + ^Ca(3)^ 

- ^ai)a3)a5) + ^a3)^a5) + '-^u^r 

32864 3752 2 5696 3 13097 

^Ca(l)Ca(5) g-Ca(3)Ca(5) - ^^Ca(5) 90~^«^"^ 

+ ^Ca(l)C(7) - ^C(3)Ca(7) + ^Ca(l)Ca(3)C(7) 
+ ^Ca(3)\„(7) - ^Ca(5)Ca(7) + ^Ca(3)C.(5)C.(7) 

95746, ,„,2 139552, , ,„,2 4921 , , , 

+ -^Ca{7r - ^^Ca(l)Ca(7)^ + ^Ca(9) 

- 2016Ca(l)Ca(9) + ^^^Ca(3)C(9) - ^Ca(l)Ca(3)Ca(9) 

75 15 
3Q872 5471 2 7Q744 

- ^Ca(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)a(9) 

= -0.008335007472438759 • • • (B.25) 
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1 11 1Q7Q 833 4Qiq 

(SfSIS^.S^,) = —- -C„(l) + ^Ca(3) - ^Ca(l)Ca(3) - ^C(3)^ 

13612,,^,^ 75923,,^, 10660,,,,,,^, 10547,,^,,,^, 

+ ^ai)a3)a5) - ^a3)^a5) - ^^5)^ 

200018^ 28217^ 290444^ ,^ _ 

- ^^Ca(l)Ca(7) + ^-Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

88312 , 59927 67088 

4^^«(2) + ^^Ca(5)Ca(7) ^Ca(3)Ca(5)Ca(7) 

162484^,^,2 234808 _^__2 3892 

" + ^^Ca(l)Ca(7)' - ^Ca(9) 

+ 3318Ca(l)Ca(9) - ^^Ca(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 

67088, ,^,2. /^N 92848, 134176 

+ ^^^«(3)\„(9) + ^^Ca(5)Ca(9) - ^^Ca(l)Ca(5)Ca(9) 

= 0.0008303788046606665 • • • (B.26) 



{S^,SIS^,SI) = ^ - gCa(l) + ^Ca(3) - ^Ca(l)Ca(3) " ^Ca(3)^ 

13612 156821 3596 5161 

~ ^35~^"^^ ~ 1800 ^"^^ Ca(l)Ca(5) + ^^Ca(3)Ca(5) 

+ ^C(l)Ca(3)C(5) - ^Ca(3)^C(5) - ^C(5)^ 
165968, ,^,2 6308, ,^,2 9584, ,^,o 10409,, 

+ ^^Ca(l)Ca(5)' + -g-Ca(3)Ca(5)' + ^Ca(5)' + ^^Ca(7) 

200543^,,,^,^, 28238,,^,^,^, 290444,,^,^,^,^,^, 

- -^Ca{l)Ca{7) + ^-Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

88312, ,^,2. 59927, 67088 , , , 

- ^^Ca(3)\„(7) + ^^Ca(5)Ca(7) - ^^Ca(3)Ca(5)Ca(7) 

162484, ,„,2 234808, , ,„,2 3892, , , 

~ + ^^Ca(l)Ca(7)' - ^Ca(9) 

+ 3318Ca(l)Ca(9) - ^^Ca(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 

75 15 
67088, , ,o, , , 92848, , , , 134176, ,,,,,,,, 

+ ^^Ca(3)\a(9) + ^^Ca(5)Ca(9) - ^^Ca(l)Ca(5)Ca(9) 

= -0.001340720075108033 • • • (B.27) 
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1 11 2S S82 7Qfil 

{S^S^.S^.S^SlSy) = ^ + — - -C(3) + ^Ca(l)Ca(3) + -J^Cai^r 

148, 90583,,^, 5887,,,,,,^, 4061,,^,,,^, 

+ -^Ca(3)^ + -^Ca(5) - ^Ca(l)Ca(5) - ^Ca(3)Ca(5) 

296^ 6716^ ,^,2^ 430573^ ,^,2 

- ^Ca(l)Ca(3)Ca(5) + ^Ca(3)'Ca(5) + -Q^Cai^f 

- '-^CMU^f - f U3)C(5)^ - '^U'^f - '-^U7) 

69223^ 13816K 47012^ 

+ -QQ-Ul)Ca{7) - -gQ-Ca(3)Ca(7) + ^^Ca(l)Ca(3)Ca(7) 

9"^? 39919 1296 

+ -fCama{7) - ^Ca(5)Ca(7) + ^Ca(3)Ca(5)Ca(7) 

34811 2 4221 ,_,,2 6543 

- ^Ca(l)Ca(9) + ^Ca(3)Ca(9) - 952C(l)Ca(3)Ca(9) 

- ^Ca(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)Ca(9) 

= -0.005996024922536831 • • • (B.28) 



{s^.sisisisisD = ^ - ^c(i) - ^C(3) + ^Ca(l)C(3) + ^Ca(3)^ 
+ + - ^C.(1)C(5) - |^C„(3)C„(5) 

- ^Ca(l)Ca(3)Ca(5) + ^C(3)^C(5) + ^C(5)^ 

- ^^Ca(l)Ca(5) - 34Ca(3)Ca(5) - -^Ca(5) - ^^Ca(7) 

11753 _ 162149 14126 _ 

+ ^^Ca(l)Ca(7) - ^^Ca(3)Ca(7) + ^^Ca(l)Ca(3)Ca(7) 

476 24209 1084 

+ ^Ca(3)^Ca(7) - ^Ca(5)Ca(7) + ^Ca(3)C.(5)Ca(7) 

+ ^^7^ - ^ai)a7)^ + ^c(9) 

- ^Ca(l)Ca(9) + ^Ca(3)C(9) - 1)C(3)C (9) 

- ^Ca(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)C(5)Ca(9) 

= 0.001715129839332883 • • • (B.29) 
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{S^.S^.S^SISISD = ^ - ^Ca(l) - ^Ca(3) + |Ca(l)Ca(3) + ^Ca(3)^ 

223,, 82489,,^, 5309,,,,,,^, 2843,,^,,,^, 

+ 1^^"^^^ + 16800^"^^^ " ^Ca(l)Ca(5) - ^Ca(3)C.(5) 

- ^Ca(l)Ca(3)Ca(5) + ^C.(3)^C.(5) + ^Ca(5)^ 

- ^C(1)C(5)^ - 34C(3)C(5)^ - ^C(5)^ - ^Ca(7) 
45563^ 160259^ 14126 ^ ,^ _ _ 

+ ^gQ-Ca(l)Ca(7) - ^gQ-Ca(3)Ca(7) + ^^Ca(l)Ca(3)Ca(7) 
+ ^Ca(3)\.(7) - ^Ca(5)Ca(7) + ^Ca(3)Ca(5)Ca(7) 
+ ^a7f - ^Ca(l)Ca(7r + 1^0(9) 

- ^Ca(l)Ca(9) + ^Ca(3)Ca(9) " ^Ca(l)Ca(3)Ca(9) 

3252 ^ /„so ^ /„s 1389 . . 2168 . , ^ ^ /„s 

- ^Ca(3)\a(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)Ca(9) 

= -0.002589603677721155 • • • (B.30) 



1 11 1537 491 317 

{S'tS-,SlS^,SISy) = ^ - ^Ca(l) + Y44^Ca(3) - — Ca(l)Ca(3) - — C„(3)^ 

- - + ^c«(i)c(5) + ^c(3)c(5) 

+ ^C(l)Ca(3)Ca(5) - ^C(3)^Ca(5) - ^Ca(5)^ 

25616, ,,,, ,^,2 ,^o. /ox. /.x2 3452,, 15553 , 

+ ^^Ca(l)Ca(5)' + 108C„(3)Ca(5)^ + -^Cai^f + ^^Ca(7) 

- ^Ca(l)Ca(7) + ^Ca(3)Ca(7) - ^C(l)Ca(3)Ca(7) 

- ^Ca(3)^Ca(7) + ^Ca(5)Ca(7) - ^C.(3)Ca(5)Ca(7) 

32809, ,„,2 12082, , 1593, , , 

- -^Ury + -^Ca{l)Ca{7f - ^Ca(9) 

+ ^C(1)C(9) - ^C(3)C(9) + ^C(1)C(3)C(9) 

+ ^C«(3)\„(9) + ^C.(5)C(9) - ^C(1)C(5)C(9) 

= -0.002089292530622660 • • • (B.31) 
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- ^C.(3)' - ^C(6) + ^C.(1)C(5) + ^C(3)C.(5) 

+^ai)c„(3)c.(5)-^c.(3)^a5)-i2f?a5)^ 

19161U,^__ 19438U 13279^ _ 

~ n^^"^^^"^^^ + -90-Ca(3)Ca(7) - -g-Ca(l)Ca(3)Ca(7) 

20174 2 2798 3068 

4^^»(3) Ca(7) + -^Ca(5)Ca(7) ^Ca(3)Ca(5)Ca(7) 

14777 _^,2 10738 _^__2 1763 

- + -^Ca(l)Ca(7)' - ^Ca(9) 

4029 12269 20174 

+ ^Ca(l)Ca(9) - ^Ca(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 

5 6 15 

3068, ,^,2. 4222,, 6136 

+ ^Ca(3)\a(9) + ^Ca(5)Ca(9) - ^Ca(l)Ca(5)Ca(9) 

= 0.0008959339234083256 • • • (B.32) 



1 13 619 229 14687 

{S^,S'^SIS:SIS^,) = ^ - ^Ca(l) + 4^C(3) - ^Ca(l)Ca(3) " ^Ca(3)^ 

79 74279 2723 5197 

- f C.(3)^ - ^C(5) + ^C(l)Ca(5) + ^C(3)C.(5) 

+ ^C(l)Ca(3)Ca(5) - ^Ca(3)\„(5) - ^Ca(5)^ 

+ ^Ca(l)Ca(5)^ + 108C(3)C„(5)^ + ^Ca(5)^ + ^Ca(7) 

- ^Ca(l)C(7) + ^C(3)C„(7) - ^C«(l)Ca(3)C(7) 

- ^Ca(3)^Ca(7) + ^Ca(5)Ca(7) - ^C.(3)Ca(5)a(7) 

32809, ,„,2 12082, , ,„,2 1593, , , 

- -^Ca(7)' + ^^Ca(l)Ca(7)' - ^Ca(9) 

+ ^C(1)C(9) - ^C(3)C(9) + ^C(1)C(3)C(9) 

+ ^C.(3)\„(9) + ^C.(5)C(9) - ^C(1)C«(5)C(9) 

= 0.001508610170462532 • • • (B.33) 
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11 2'i 1 fi 121 

(S^SySISlS^Sy) = -C(l) + — Ca(3) - ^Ca(l)Ca(3) - ^Ca(3)^ 

73 , ,^,o 1541 413,,,,,,^, 539,,^,,,^, 

- T35^«(^) + 50400^«(^) + ^C(1)C(5) - ^C(3)C„(5) 

+ ^Ca(l)Ca(3)Ca(5) - |^Ca(3)^Ca(5) - ^Ca(5)^ 

+ ^Ca(l)C(5)^ + f Ca(3)C(5)^ + ^C(5)^ - ^C(7) 
217 2578 1477 

- ^Ca(l)Ca(7) + ^Ca(3)Ca(7) - ^Ca(l)Ca(3)Ca(7) 

- ^Ca(3)\.(7) + ^Ca(5)Ca(7) - ^Ca(3)Ca(5)Ca(7) 

- ^a7f + ^ai)Ca(7f + ^Ca(9) 

1 flS 801 7 448 

+ ^Ca(l)Ca(9) - ^Ca(3)Ca(9) + — Ca(l)Ca(3)Ca(9) 
+ ^Ca(3)\„(9) + ^Ca(5)Ca(9) - ^Ca(l)Ca(5)Ca(9) 

= 0.0007223711440396728 • • • (B.34) 



1 3'i3 5fiQ 327 

(S'^S'^SlStSlS^,) = ^ + -Ca(l) - ^Ca(3) + ^Ca(l)Ca(3) + — Ca(3)^ 

ai^r + |^Ca(5) - ^Ca(l)Ca(5) " ^Ca(3)Ca(5) 

a(l)C(3)C(5) + ^C(3)\.(5) + ^Ca(5)^ 

Ca(l)Ca(5)^ - ^Ca(3)Ca(5)^ - ^C(5)^ - ^Ca(7) 

22497 341 fi 

Ca(l)Ca(7) - ^Ca(3)Ca(7) + ^Ca(l)Ca(3)Ca(7) 

Ca(3)^Ca(7) - ^Ca(5)Ca(7) + ^Ca(3)Ca(5)Ca(7) 

C(7)^-^C(l)a7)^ + ^C(9) 
1 5Qfil IfT^fi 

a(l)Ca(9) + ^Ca(3)Ca(9) - ^Ca(l)C(3)Ca(9) 
a(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)Ca(9) 

0.001419281470187951 • • • (B.35) 
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(S^.S^.SlStSISl) = ^ + J^Ca(l) - ^C(3) + ^Ca(l)Ca(3) + ^Ca(3)^ 
+ ^Ca(3)^ + ^C(5) - ^Ca(l)Ca(5) - ^C(3)C(5) 

- ^Ca(l)Ca(3)a(5) + ^Ca(3)^Ca(5) + ^Ca(5)^ 

- ^C(1)C(5)^ - ^C(3)C(5)^ - ^C(5)^ - ^C(7) 

+ ^Ca(l)Ca(7) - ^Ca(3)Ca(7) + ^Ca(l)Ca(3)Ca(7) 

+ ^Ca(3)\.(7) - ^Ca(5)Ca(7) + ^Ca(3)Ca(5)Ca(7) 

7609 2 2758 2 251 

+ -^U7f - -^Ca(l)Ca(7)' + ^Ca(9) 

- ^Ca(l)Ca(9) + ^Ca(3)Ca(9) - ^Ca(l)Ca(3)Ca(9) 

- ^Ca(3)\„(9) - ^Ca(5)Ca(9) + ^Ca(l)Ca(5)Ca(9) 

= -0.002041037236271903 • • • (B.36) 



1 7 116Q Q41 3391 

1474 3 150971 5662 1237^ , , 

- 135^^^'^ - 16800 ^^^'^ + ^C.(1)C(5) + ^C(3)C(5) 

+ ^Ca(l)Ca(3)Ca(5) - ^Ca(3)\a(5) - ^Ca(5)^ 

5896, ,,,, ,^,2 671^,^,^,^,2 7148,, 21893,, 

+ ^Ca(l)Ca(5)' + -g-Ca(3)Ca(5)' + ^Ca(5)' + ^^Ca(7) 

22673,,,,,,^, 5013,,^,,,^, 10318,,^,,,^,,,^, 

- ^^Ca(l)Ca(7) + ^Ca(3)Ca(7) - ^^Ca(l)Ca(3)Ca(7) 

9394 2179 7148 

- ^C.(3)^Ca(7) + ^Ca(5)Ca(7) - ^Ca(3)Ca(5)Ca(7) 

17059, ,„,2 25018, , ,„,2 4441, , , 

1917 71282 9394 

+ ^Ca(l)Ca(9) - ^C(3)C(9) + ^C(1)C(3)C(9) 

7148 9748 14296 

+ ^Ca(3)\.(9) + ^Ca(5)Ca(9) - ^Ca(l)Ca(5)Ca(9) 

= -0.0001528594146694370 • • • (B.37) 
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529, 134587,,^, 6328,,,,,,^, 2737,,^,,,^, 

" 1^^'^^^^ " 12600"^"^^^ + ^Ca(l)Ca(5) + ^Ca(3)Ca(5) 
+ ^Ca(l)Ca(3)C.(5) - ^Ca(3)^Ca(5) - ^Ca(5)^ 

+ ^C(1)C(5)^ + 82C(3)C(5)^ + |^Ca(5)^ + ^C(7) 

- ^Ca(l)Ca(7) + ^Ca(3)Ca(7) - ^Ca(l)Ca(3)Ca(7) 

- ^Ca(3)\a(7) + ^Ca(5)Ca(7) " ^Ca(3)Ca(5)Ca(7) 

2086, _2 3052,, 2 4581 

- -^Ca{7f + -^Ca{l)Ca{7f - ^Ca(9) 

4121 26024 3444 

+ ^Ca(l)Ca(9) - ^Ca(3)Ca(9) + ^Ca(l)Ca(3)Ca(9) 
+ ^Ca(3)\„(9) + ^Ca(5)Ca(9) - ^Ca(l)Ca(5)Ca(9) 

= -0.0003767104281403536 • • • (B.38) 
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